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Abstract 

We review the work done so far aimed at modeling in an alternative way the dark mat- 
ter in the Universe: the scalar field/ Bose-Einstein condensate dark matter (SFDM/BEC) 
model. We discuss a number of important achievements and characteristics of the model. 
We also describe some of our most recent results and predictions of the model compared 
to those of the standard model of ACDM. 

1 Introduction 

It is a pleasure for us to review the different theoretical basis of the SFDM/BEC model as 
a dark matter (DM) candidate. We think this review is important for mainly two reasons: 
the considerable progress in the model since it was proposed as a serious candidate to the 
dark matter paradigm, improved theoretical understanding of the nature of the DM and 
the significant advances in tlic cosmological and astronomical observations are leading us 
to put more constraints and will allow us to test the model and decide if it can still stand 
as a viable DM paradigm or if it should be discarded. 

Recent observations of the Universe have found that only 4 percent of the total content 
of the Universe is baryonic matter, being 22 percent of the rest remaining non-baryonic 
dark matter (DM) and the rest in some form of cosmological constant. 

The incorporation of a new kind of DM different from that proposed by tlic standard 
model, also known by A-Cold Dark Matter model (ACDM), into the Big Bang Theory 
holds out the possibility of giving alternative answers to some of the unsolved issues of 
the standard cosmological model. Several authors have proposed interesting alternatives 
in where they try to solve the difficulties that the ACDM scenario seems not to solve. In 



fact, the alternative Scalar Field Dark Matter (SFDM/BEC) scenario has received much 
attention in the last few years. The main idea is simple, in this models the nature of the 
DM is completely determined by a fundamental scalar field $ [2] . 

The idea was first considered in [5H] and independently introduced in [331 1351 [53] sug- 
gesting bosonic dark matter as a model for galactic halos, see also [35]. In the SFDM/BEC 
model, DM halos can be described, in the non-relativistic regime, as Newtonian gravita- 
tional condensates made up of ultra-light bosons that condense into a single macroscopic 
wave function. 

Several authors have introduced a dynamical scalar field with a certain potential 
as a candidate to be the dark matter, although there is not yet an agreement for 
the correct form of the potential field. Other interesting works consider a single scalar 
field to unify the description of dark matter, dark energy and inflation [70 l l7 H l54]. 

Different issues of the cosmological behavior of the SFDM/BEC model have been 
studied in a wide variety of approaches, see for example [HUl [21 US] [SDl 11181 [751 [771 11061 
[TU71 HH [Sni [HIMl [SlIMl [Ml [US]- For example, [5^ proposed fuzzy dark matter 
composed of ultra-light scalar particles initially in the form of a BEC. Recently [501 [52] 
developed a further analysis of the cosmological dynamics of SFDM/BEC as well as the 
evolution of their fluctuations (see also [26]). In the same direction, [ll4l [76] studied the 
growth of scalar fluctuations and the formations of large-scale structure within a fluid 
and a field approach for the SFDM/BEC model. 

In addition, many numerical simulations have been performed to study the gravita- 
tional collapse of the SFDM/BEC model [HHHHaililiaiiaEHlIHHlIHSIISl- [HHS] 
found an approximate analytical expression and numerical solutions of the mass-radius 
relation of SFDM/BEC halos. Recently, [103] gave constraints on the boson mass to 
form and maintain more than one vortex in SFDM/BEC halos. These constraints are 
in agreement with the ultra-light mass found in previous works (see for example, |59|). 
Lately, [TSj performed N-body simulations to study the dynamics observed in the Ursa 
Minor dwarf galaxy. They modeled the dark matter halo of Ursa Minor as a SFDM/BEC 
halo to establish constraints for the bosons mass. Moreover, they introduced a dynamical 
friction analysis within the SFDM/BEC model to study the wide distribution of globular 
clusters in Fornax. An overall good agreement is found for the ultra-light mass of bosonic 
dark matter. 

In this paper we review a number or results and important characteristics of the 
SFDM/BEC model, its dynamical mechanism and some of its predictions. We also discuss 
some of the trending topics nowadays in the subject which attempt to predict and ask 
how well the model is achieving its goals. As we will see, a number of them studies results 
which are in reasonable agreement with the general features required by the theory and 
the data. 

The outline of the paper is as follows. In Section 1 we have given a brief synopsis 
on the current state-of-the-art of the model. In the Section 2 we describe the dark 
matter paradigm and briefiy resume the standard model of cosmology. In Section 3 and 
4 we describe in some detail why Scalar Field/Bose-Einstein condensate Dark Matter 
(SFDM/BEC) can be a good alternative candidate to be the dark matter in our Universe, 
we summarize some representative papers for this sections, and in Section 5 we include 
topics for future works and our conclusions. 
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2 The Dark Matter Paradigm and the Standard Model 
of Cosmology 

In this millennium, new technologies are opening wider windows to explore our Universe. 
For some time we could only relay on inaccurate evidence found in the local neighborhood 
of our galaxy to infer the history of our Universe, now it turns it is possible for us to 
see the evolution of the Universe as far as 100,000 years after the Big Bang and in more 
detail. 

With these advances, nowadays some inquires of our cosmic evolution can be deter- 
mined by giving an answer to question like: How much matter is in the Universe? 

Since the discovery of the expansion of the Universe done by Hubble and Slipher [5B] 
in the 1920's, the common believe had been that all energy in the Universe was in the 
form of radiation and ordinary matter (electrons, protons, neutrons, etc.). Over the past 
few decades, theories concerning the stability of galaxies ( |108| [9]) indicated that most 
of the mass in our Universe is dark (i.e., it does not emits or absorbs light [SOI 
therefore resulting unobservable by telescopes. The suggestion that "dark matter" may 
form a large fraction of the density in the Universe was raised by Zwicky in 1937. Back 
then he used the virial theorem to obtain the average mass of galaxies within the Coma 
cluster and obtained a value much larger than the mass of the luminous material, he then 
realized that some mass was "missing" in order to account for observations. This missing 
mass problem was confirmed many years later by more accurate measurements of rotation 
curves of disc galaxies, ^191 [Ml [TU51 The rotation curves of neutral hydrogen clouds 
in spiral galaxies measured by the Doppler effect are now found to be roughly flat with 
a typical rotation velocity equal to Voo ^ 200km/s up to the maximum observed radius 
of about 50kpc. With these observations the mass profile results much more extended 
than the common distributions which typically converge within ^ lOkpc. This would 
imply that galaxies might be surrounded by an extended halo of dark matter whose 
mass M(r) rv'^/G increases linearly with radius (here r is the radius and G Newton's 
gravitational constant) . 

In the 1980's, the proposal of dark matter found its basis in the so called "inflationary 
scenario" [HJ [75J [T], a theory of the first 10^"^° seconds developed to give answer to 
several questions left unanswered by the Big Bang model, like for example: Why is the 
Universe so homogeneous and isotropic? and; Where did the initial homogeneities that 
gave rise to the structures we see today came from? |43l [121 Hll 11131 [53] . The inflationary 
theory predicts that the Universe is spatially flat; which according to Einstein's theory 
of general relativity, this fixes the total energy density of the Universe making it equal to 
the critical value, pc = 3Hq/8ttG ^ 1.7 x 10~^^g cm^, where Hq is the current value of 
the Hubble parameter. 

Several astrophysical observations of distant type la supernovae have also shown that 
the content of the Universe is made of about 74% of dark energy, 22% of dark matter 
and 4% of baryonic (visible) matter, [HTl [H [TTTl [TUni [SZl [H] . Inflation thus seemed to 
call for dark matter. 

It then results that most matter in the Universe is non-luminous. The observed flat- 
ness of the galactic rotation curved indicating the presence of dark matter halos around 
galaxies IIOSJ . Summed to these evidences are the observations of the cosmic microwave 
background (CMB) anisotropics ( (112t l57j) combined with large-scale structure and type 
la supernova luminosity data ( |100l [TUTl 11021 which all together constrain the cos- 
mological parameters also finding once more that visible matter contributes only about 
4% of the energy density of the Universe, gravitational lensing [M] and X-ray spectra 
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( [69l [37l 11161 [36l [87]) in elliptical galaxies, and the high velocity dispersion and gas 
temperature in clusters of galaxies ( |116[ [74]). all of them leading to a picture in which 
galaxies are composed of a luminous galactic disk surrounded by a galactic halo of dark 
matter. Also the relative contribution of the dark matter component is usually specified 
in terms of the mass-to-light ratio, M/L; which reflects the total amount of mass relative 
to the total light within a given scale. The increase on this ratio suggests that there is 
relatively more dark than luminous matter with increasing scale [108j . This has led to 
the general belief that clusters have more dark matter per unit luminosity than individual 
galaxies and that superclusters may have even more. This widely accepted monotonic 
increase of M/L with scale determines to a large extent the prevalent views about the 
location of the dark matter and the total mass density of the Universe. Recent studies of 
the dependence of the mass-to-light ratio on scale indicate that M/L is nearly constant 
on large scales ranging up to supercluster size (lOMpc), suggesting no additional dark 
matter is tucked away on large scales [TD]. More recently, a clear separation between the 
center of baryonic matter and the total center of mass was observed in the Bullet cluster 
([S7]) and later in other galaxy cluster collisions ([5U]). 

The observational evidence for dark matter continues to grow, and particle physicists 
have proposed various particles, motivated by supersymmetry and unified theories, that 
could reasonably explain it. These observations reinforce the claim that dark matter 
is indeed composed of weakly interacting particles and is not a modification of gravity. 
However, even taking into account all these results the properties of dark matter are still 
mysterious. 

It then also results that an important question in cosmology has to do with knowing 
the nature of the so far undefined one quarter part of the content of the Universe, the 
dark matter. As mentioned before, the cosmological observations seem to support the 
idea that dark matter can be made of some kind of non-baryonic, non-relativistic and 
weakly interactive massive particle. Many efforts trying to give an answer to this question 
have been done in the past few decades, mainly motivated by the idea that the answer 
will probably change our understanding of the Universe and its dynamics. One of the 
explanations for DM is the SFDM model. 

In the Standard Model of cosmology, the total energy density of the Universe is domi- 
nated today by the densities of two components: the "dark matter" which has an attrac- 
tive gravitational effect like usual matter and the " dark energy" which can be considered 
as a kind of vacuum energy with a negative pressure, which seems constant today (i.e., 
a cosmological constant. A). Although the real nature of these two components remains 
unknown, in the standard model dark matter is generally modeled as a system of collision- 
less particles. This is known as the "A Cold Dark Matter" model, which predicts that the 
Universe contains primarily cold neutral weakly interactive massive particles (WIMPs) 
which are non-baryonic (|92 |I123) ). pressureless, behave like a cold gas, one beyond those 
existing in the Standard Model of particle physics and have not yet been detected in 
accelerators or specialized indirect searches, in particular, the lightest supersymmetric 
particles, the most popular of which is the neutralino, with a particle mass of the order 
of lOOGeV. Efforts are underway to measure the presence of these particles, but no direct 
detection has yet been reported. 

In order to explain observational data, the ACDM model was developed, [MIISI] and it 
is the most simple possibility. In the Standard Model of cosmology the matter component 
f^Af ^ 26% of the Universe decomposes itself into baryons, neutrinos, etc., and cold 
dark matter which would be responsible for the formation of structure in the Universe. 
Observations indicate that stars and dust (baryons) represent something like 0.4% of the 
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total content of matter in the Universe. The measurements of neutrino masses indicate 
that these contribute nearly with the same amount as matter. In other words, flM = 
f2m + ^DM = flb + ^u + ■■■ + ^DM ~ + ^CDAU where i^cDM represents the cold dark 
matter part of the matter contributions, and has a value of ficD a/ ^ 0.22. The value of the 
amount of baryonic matter is in accordance with the limits imposed by nucleosynthesis, 
[llOj . This model then considers a flat Universe (J7a + Hm = 1) with 96% of unknown 
matter but which is of great importance in the cosmological context. It also supposes 
a homogeneous and isotropic Universe which evolution can be best described today by 
Friedmann's equations coming from general relativity and whose main ingredients can be 
described by fluids with characteristics very similar to those we see in our Universe. We 
now know that the Universe is not exactly homogeneous and isotropic, but the standard 
model does give a framework within which the evolution of structures such as galaxies 
or clusters of galaxies can be studied with their origins coming from small fluctuations 
in the density of the early Universe. The model assumes a " scale- invariant" spectrum 
of initial density fluctuations, a spectrum in which the magnitude of the inhomogeneity 
is the same on all length scales, again as predicted by standard inflationary cosmology 
[l2l l42l 11131 1121j . Moreover, ACDM seems to be until today the most successful model 
fitting current cosmological observations fVD . 

The ACDM model successfully describes the accelerated expansion of the Universe, 
it explains the Cosmic Microwave Background radiation in great detail and provides a 
framework within which one can understand the large-scale isotropy of the Universe, 
it also describes the important characteristics of the origin, nature and evolution of the 
density fluctuations which are believed to give rise to galaxies and other cosmic structures, 
the Lyman-a forest, the large scale matter distribution, and the main aspects of the 
formation an the evolution of virialized cosmological objects. So far the ACDM model is 
consistent with the observed cluster abundance at z ~ 0, it then predicts a relatively little 
change in the number density of rich clusters as a function of redshift because, due to the 
low matter density, hardly any structure growth has occurred since z ^ 1 . The ACDM 
model can then be "forced" to agree approximately with both the cluster abundance on 
small scales and the CMB fluctuations on large scales by tilting the power spectrum (by 
about 30%) from its standard shape. This tilted variant of the ACDM model then results 
nearly consistent with observations. The power spectra of ACDM can then be normalized 
so that it agrees with both the CMB and cluster observations. But as the estimates of 
the cold dark matter density become more precise, it becomes even more imperative for 
its composition to be identified. 

There remain, however, certain conflicts at galactic scales, like the cusp profile of 
central densities in galactic halos, the overpopulation of substructures predicted by N- 
body numerical simulations which are an order of magnitude larger than what has been 
observed, among others, see for example [27 t [6T | \8E \ I95j. And until today the nature of 
the dark matter that binds galaxies remains an open question. 

3 Why Scalar Field Dark Matter? 

In the big bang model, gravity plays an essential role: it collects the dark matter in 
concentrated regions called 'dark matter halos'. Within these large dark matter halos, 
the baryons are believed to be so dense that they radiate enough energy to collapse into 
galaxies and stars. The most massive halos, hosts for the brightest galaxies, are formed in 
regions with the highest local mass density. Less massive halos, hosts for the less bright 
galaxies, appear in regions with low local densities f94]. These situations appear to be 
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the same as in our extragalactic neighborhood, but there are still problems. Despite all 
its successful achievements the ACDM model requires further considerations. 

The ACDM paradigm faces several challenges to explain observations at galactic 
scales, such as the central densities of dark halos, dwarfs and Low Surface Brightness 
(LSB) galaxies, the excess of satellite galaxies predicted by N-body simulations, the for- 
mation of bars in disc galaxies, etc.[HHl HZl ES]- In other words, there is not a match 
between ACDM predictions at galactic scales and what is being observed. Problems with 
an otherwise successful model are often the key to a new and deeper understanding. 

Observations point out to a better understanding of the theory beginning with the 
Local Void, which contains just a few galaxies that are larger than expected. This problem 
would be solved if structure grew faster than in the standard theory, therefore filling the 
local void and giving rise to more matter in the surroundings (|94j). 

Another problem arises for the so called pure disk galaxies, which do not appear in 
numerical simulations of structure formation in the Standard Model. These problems 
would be solved again if the structure grew faster than it does in the standard paradigm 

m- 

On the other hand, [66] also found that the collision velocity of 3000km/s at R200 
for the Bullet Cluster is very unlikely within the ACDM paradigm, which moves it to a 
challenge for the Standard Model of cosmology. 

A final example of inconsistencies can be seen in a paper by [115] , who found anomalies 
in the mass power spectrum obtained by the SDSS and the one obtained with the ACDM 
model, i.e., anomalies in the predicted large-scale structure of the Universe. With these 
and other results it seems necessary to change the ACDM paradigm to try and explain 
the formation of structure in the Universe. 

Given these discrepancies, it seems necessary to explore alternatives to the paradigm 
of structure formation. These are some of the reasons why we need to look for alternative 
candidates that can explain the structure formation at cosmological level, the observed 
amount of dwarf galaxies, and the dark matter density profiles in the core of galaxies. 
Recently, several alternative models have been proposed. 

One of them invokes a scalar field as dark matter in the Universe [5TJ [52]. This 
model supposes that dark matter is a real or complex scalar field <& minimally coupled 
to gravity, endowed with a scalar potential V{^) and that at some temperature it only 
interacts gravitationally with the rest of the matter. This scalar field can be added to the 
particles standard model lagrangian or to the general relativity one, supposing that the 
coupling constant with the rest of the matter is very small. It has been also suggested 
that this scalar field can be derived from higher dimensional theories. It has also been 
proposed that this dark matter scalar field, i.e., this spin-0 fundamental interaction, 
could lead to the formation of Bose-Einstein condensate in the way of cosmic structure 
( [8OI [33l [50I [52I [55] ) with an ultra-light mass of order m ~ 10~^^eV. From this mass 
it follows that the critical temperature of condensation T^, ^ l/m^/'^'^TeV is very high, 
therefore, they may form Bose-Einstein condensate drops very early in the Universe [ST] 
that behave as cold DM. Lee and Koh [53], and independently Matos and Guzman [ii] . 
suggested bosonic dark matter as a model for galactic halos. In addition, the Compton 
length Ac = 2'Kh/ra associated to this boson results of about ~kpc, and corresponds to 
the dark halo size of typical galaxies in the Universe. Thus, it has been sopported that 
these drops are the halos of galaxies (see [HO]), i-e., that halos are huge drops of SF. In 
a recent paper, Urefia jll8| studied the conditions for the formation of a SFDM/BEC in 
the Universe, also concluding that SFDM/BEC particles must be ultra-light bosons. 

In the SFDM model the initial halos of galaxies do not form hierarchically, they are 
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formed at the same time and in the same way when the Universe reaches the critical 
temperature of condensation of the SF. From this it follows that galaxies can share some 
properties because they formed in the same manner and at the same moment |67j . There- 
fore, from this paradigm we have to expect that there exists well formed galaxy halos at 
higher redshifts than in the ACDM model. Recently Suarez and Matos |114j developed 
a hydrodynamical approach for the structure formation in the Universe with the scalar 
potential V{^) = m^<i>^/2 + A$'^/4. They found that when A = the evolution of per- 
turbations of the SFDM model compared to those of ACDM are identical. They also 
showed that this potential can lead to the early formation of gravitational structures in 
the Universe depending on the sign of the self-interaction parameter A. 

The most simple model having both an exponential behavior and a minimum is a 
cosh-like potential. Another interesting work was done by |78j an independently by |109) 
who used a potential of the form V{^) = Vb[cosh(?7$ — 1)] where Vq and r] are constants 
to explain the core density problem for disc galaxy halos in the ACDM model (see also 
[551 and to perform the fist cosmological analysis in the context of SFDM. They 
showed that the evolution of the Universe, its expansion rate and the growth of linear 
perturbations in this model are identical as those derived in the standard model. In 
|44j they model the dark matter in spiral galaxies, assuming dark matter as an arbitrary 
scalar field endowed with a scalar potential. 

Another scalar potential widely used to describe dark matter is V{^) — |1171 
I5T] . This potential is very interesting because it can mimic the cosmological evolution of 
the Universe predicted by the ACDM model. Also, it is known that an exponential-like 
scalar field potential fits very well the cosmological constraints due to the form of its 
solutions (see for example [^112211^ '). If the self-interaction of the SF is considered, we 
need to add a quartic term to the SF potential [31 HI 13 151 HI] , in this case the equation 
of state of the SF results to be that of a polytrope of index n=l (see [1141 [50115^ IHIS5]). 

Another interesting result is that the predicted density of neutrinos at the recombina- 
tion epoch is in agreement with the observations of the Wilkinson Microwave Anisotropy 
Probe (WMAP). In the same direction, Rodn'guez-Montoya et al. |106j studied ultra- 
light bosons as dark matter in the Universe with the framework of kinetic theory, through 
the Boltzmann-Einstein equations, and they found that this kind of ultra-light particles 
is consistent with the acoustic peaks of the cosmic microwave background radiation if the 
boson mass is around m ~ 10~^^eV. 

[64] pointed out that SFDM/BEC can explain the spatial separation of the dark 
matter from visible matter, as derived from X-ray maps and weak gravitational lensing, 
in the Bullet Cluster, see also [?7] . 

Other works have used the bosonic dark matter model to explain the structure for- 
mation via high-resolution simulations. [1181 1119) reviewed the key properties that may 
arise from the bosonic nature of SFDM models. On the other hand, several authors have 
numerically studied the formation, collapse and viralization of SFDM/BEC halos as well 
as the dynamics of the SFDM around black holes [H iZlllSl [HI [13 [Ml 133 (Hj ■ In [2] , 
Alcubierre et al. found that the critical mass for collapse of the SF is of the order of a 
Milky Way-sized halo mass. This suggests that SFDM/BEC can be plausible candidate 
to dark matter in galactic halos. In addition, Lora et al.,[73], studied, through N-body 
simulations, the dynamics of Ursa Minor dwarf galaxy and its stellar clump assuming a 
SFDM/BEC halo to establish constraints for the boson mass. Moreover, they introduced 
a dynamical friction analysis with the SFDM/BEC model to study the distribution of 
globular clusters in Fornax. An overall good agreement is found for the ultra-light mass 
^ 10~22gY bosonic dark matter. 
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In this model the scalar particles with ultra-light mass are such that their wave prop- 
erties avoid the cusp problem and reduce the high number of small satellites by the 
quantum uncertainty principle ( [551 [751 [511 [CTl [TM] ). Robles and Matos [104] . (see also 
[T51 [TBI [5T] ) showed that BEC dark matter halos fit very well high-resolution rotation 
curves of LSB galaxies, and that the constant density core in dark halos can be repro- 
duced. Also, [ST] showed how the SFDM/BEC paradi gm is a good alternative to explain 
the common mass of the dark halos of dwarf spheroidal galaxies. Recently, Rindler-Daller 
and Shapiro, jl03j, investigated the formation of vortex in SFDM/BEC halos. They found 
constraints on the boson mass in agreement with the ultra-light mass found in previous 
works (see also [5^1125) ). some of these issues will be discussed in more detail in Section 
4. 

Summarizing, it is remarkable that with only one free parameter, the ultra-light scalar 
field mass (m ~ IQ-^^eV), the SFDM model fits: 

i) The evolution of the cosmological densities [HI] • 

ii) The central density profile of the dark matter is flat [IB] . 

iii) The acoustic peaks of the cosmic microwave background [106j . 

iv) The scalar field has a natural cut off, thus the substructures in cluster galaxies 
are suppressed naturally. With a scalar field mass of to$ ~ 10~^^eV the amount of 
substructures is compatible with the ones observed [5D1 [551 1114] . 

v) We expect that SFDM forms galaxies earlier than the ACDM model, because they 
form BEC's at a critical temperature Tc >>MeV. So if the SFDM model is right, we have 
to see big galaxies at high redshifts with similar features [HH I114j . 

vi) Adding self-interaction and Temperature correcctions, the rotation curves of big 
galaxies and LSB galaxies [Ml [MHHiiil [S] . 

vii) With this mass, the critical mass of collapse for a real scalar field is just lO^^M©, 
i. e., the observed in galaxy halos [2]- 

vii) The observed properties of dwarf galaxies, i. e., the minimum length scale, the 
minimum mass scale, and their independence from the brightness [67j . 

ix) And recently it has been demonstrated that the SFDM halos would have cores 
large enough to explain the longevity of the cold clump in Ursa Minor and the wide 
distribution of globular clusters in Fornax [75] . 

Then, the SFDM/BEC model has provided to be a good candidate for dark matter 
halos of galaxies in the Universe because it can explain many aspects where the standard 
model of cosmology fails ([ffl HHIIMI EHIMl ESI ISH [lQ4]). Therefore, not only the 
many successful predictions of the Standard Model of cosmology at large scales are well 
reproduced by SFDM, but also the ones at galactic scales. The scalar field models presents 
some advantages over the standard ACDM model like the ones mentioned above. Also, 
its self-interaction can, in principle, explain the smoothness of the energy density profile 
in the core of galaxies [80l [84]. Nevertheless, its important to remark that when a new 
dark matter candidate is proposed the study of the final object that will be formed as a 
result of a gravitational collapse is always an important but difficult task that requires 
continues work since the baryonic physics is still poorly understood. 

4 Current Status of the SFDM model 

4.1 Self-gravitating Bose-Einstein condensate dark matter 

Following [TS], dark matter halos as a self-gravitating Bose-Einstein condensate with 
short-range interactions have been widely discussed [1051 11041 [Ml [55] . In these models. 
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it is supposed that the cosmic BEC has a relatively low mean mass density so that the 
Newtonian approximation can be used. In the literature, when T = all the bosons have 
condensed and the system can be described by one order parameter ■(/'(r, i), called the 
condensate wave function. In the mean-field approximation, the ground state properties 
of the condensate are described by the Gross-Pitaevskii equation: 

d->b 

ih^{f,t) = - — Ail;{f,t) + m(S>tot{r,t)i^{r,t). 
ot 2m 

where $tot is the total potential exerted on the condensate. With this equation, [23 stud- 
ied the structure and the stability of a self-gravitating BEC with short-range interactions. 
In this case, the results obtained in the absence of self-coupling and the results of [18] 
obtained for self-coupled BECs in the Thomas-Fermi approximation have been connected. 
The case of attractive short-range interactions where considered and the existence of a 
maximum mass above which no equilibrium state exists was found. 

This study was motivated by the proposal that dark matter halos could be gigantic 
cosmic BEC's, [55j |63l HQI [SJ [18]. In this case, gravitational collapse is prevented by 
the Heisenberg uncertainty principle or by the short-range interaction. This suggestion 
still remains highly speculative since the nature of dark matter remains unknown. On 
the other hand, whatever the nature of its constituents, if dark matter is viewed as a 
coUisionless system described by the Vlasov equation, dark matter halos could result 
from processes of violent coUisionless relaxation. In that case, gravitational collapse 
can be prevented by a Lynden-Bells type of exclusion principle, because this form of 
relaxation could result more rapid and efficient than a collisional relaxation. Furthermore, 
it generates a density profile with a flat core and a outer density profile for the 
halo, yielding fiat rotation curves. These features resulted remarkably consistent with 
observations making this alternative scenario quite attractive. 

In [25] the same author obtained the exact mass-radius relation of self-gravitating 
BECs with short-range interactions by numerically solving the equation of hydrostatic 
equilibrium this time taking into account quantum effects. He compared his results with 
the approximate analytical relation obtained in |24j from a Gaussian ansatz. He found 
that the Gaussian ansatz always provided a good qualitative agreement with the exact 
solution, and that the agreement was quantitatively very good. 

In one of his most recent works Chavanis [25] assumed that the dark matter in the 
universe could be a self-gravitating BEC with short-range interactions, and he then the- 
oretically explored the consequences of this hypothesis. This time he considered the 
possibilities of positive and negative scattering lengths. 

At the level of dark matter halos a positive scattering length, equivalent to a repulsive 
self-interaction generating a positive pressure, is able to stabilize the halos with respect 
to gravitational collapse. This leads to dark matter halos without density cusps with an 
effective equation of state equal to that of a polytrope of index n = 1. Alternatively, if 
the scattering length is negative, equivalent to an attractive self-interaction generating a 
negative pressure, the dark matter results very unstable and collapses above a very small 
critical mass M„iax ~ l.012h/ y/\a^\Grn, where Os is the scattering length. When these 
ideas where applied to an infinite homogeneous cosmic fluid, it was found that a negative 
scattering length can increase the maximum growth rate of the instability and accelerate 
the formation of structures. The virtues of these results could be combined by assuming 
that the scattering length changes sign in the course of the evolution. It could be initially 
negative to help with the formation of structures and become positive to prevent complete 
gravitational collapse. The mechanism behind this change of sign remains unknown. 
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However, some terrestrial experiments have demonstrated that certain atoms can have 
negative scattering lengths and that it is possible in principle to manipulate the value 
and the sign of . 

It has also been found that a SFDM/BEC universe with positive scattering length, 
having a positive pressure, is not qualitatively very different from a classical Einstein-de 
Sitter universe. It also emerges at a primordial time t = from a big-bang singular- 
ity where the density results infinite, and then undergoes a decelerating expansion. A 
difference, however, is that the initial scale factor a{0) is finite. On the other hand, a 
SFDM/BEC universe with an always negative scattering length, having a negative pres- 
sure, markedly differs from previous models. It starts from t — > — oo with a vanishing 
radius and a finite density, it has an initial accelerating expansion then decelerates and 
asymptotically behaves like the EdS universe. This model universe exists for any time 
in the past and there is no big-bang singularity. When the effect of radiation, baryonic 
matter and dark energy where added, the picture resulted quite different. In that case, 
a SFDM/BEC universe with attractive or repulsive self-interaction started from a singu- 
larity at t — where the density was infinite. It first experiences a phase of decelerating 
expansion followed by a phase of accelerating expansion. For k — > {k being the wave 
number) the standard ACDM model was recovered but for k 0, the evolution of the 
scale factor in a SFDM/BEC universe turned out substantially different. The model with 
k > expanded more rapidly than the standard model. The initial scale factor is finite 
and the radiation never dominates. The model with fc < expanded less rapidly than 
the standard model. The initial scale factor vanishes and the radiation dominates leading 
to a decelerating expansion. In both models, the dark energy dominates at large times 
leading to an accelerating expansion. Finally, a dark fluid with generalized equation of 
state p — {ap + kp'^)(? having a pressure component p — kp^c^ similar to a BEC dark 
matter and a component p = apc^ mimicking the effect of the cosmological constant were 
considered. Optimal parameters (a, k) that gave a good agreement with the standard 
model where found. Also the growth of perturbations in these different models where 
studied and confirmed the previous observation of Harko, [52) . that the density contrast 
increases more rapidly in a BEC universe than in the standard model. 

In conclusion, it was pointed out that the idea that dark matter could be a BEC 
results fascinating and probably deserves further research. 

4.2 BEC dark matter and cosmological perturbations 

Through out this work we have shown how the SFDM/BEC model could be a serious 
alternative to the dark matter in the Universe. In j76) they studied in quite some detail 
the growth and virialization of $^-dark matter perturbations in the linear and nonlinear 
regimes. Following the spherical collapse model, they also studied the nonlinear regime 
of the evolution of <I>^-dark matter perturbations. They showed that the evolution of 
an overdense region of ^^-dark matter can collapse and virialize in a bound structure. 
However, they found that the scalar perturbations collapse at earlier times of the Universe 
than those in the CDM model. Thus, the standard and the SFDM/BEC model can be 
confronted in their predictions concerning the formation of the first galaxies. Massive 
galaxies at high redshifts is a prediction of the model and may be used to distinguish 
between SFDM/BEC paradigm and CDM. 

As the study in |76j is detail in the analysis of perturbations we include a brief sum- 
mary. In the study of the cosmological dynamics of the SFDM model it was considered 
the simplest case: a single scalar field $(2:, t), with self-interacting double-well potential. 
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They used the potential 



In a very early stage of the Universe, this scalar field was in local thermodynamic equi- 
librium with its surroundings see |21j . At some time, the scalar field decoupled from the 
rest of the matter and started a lonely journey with its temperature T decreased by the 
expansion of the Universe. Thus, it is considered the scalar field in a thermal bath of 
temperature T, whose scalar field potential, extended to one loop corrections, is given by 

where m is a mass parameter before the breaking of symmetry and A is the self-interacting 
constant. From here, it can be calculated the critical temperature Tc at which the Z2 
symmetry of the real SF breaks. To do that, they calculated the critical points of the 
scalar potential from 

= I -m^ + A$2 + ^T^ 1 $. 



The negative term —iv? permits the breaking of symmetry of the potential. One critical 
point is at $ = 0. If the temperature T is high enough, the scalar potential has a 
minimum at this critical point. Furthermore, the critical temperature in which $ = 
becomes a maximum is 



2 



4m 



A 

This critical temperature defines the symmetry breaking scale of the scalar field. 

To study the dynamics of the SFDM in the background Universe it is assumed a 
Friedmann-Lemaitre-Robertson- Walker metric with scale factor a{t). The background 
Universe was composed by SFDM ($0) endowed with a scalar potential V = V{^q), 
baryons (b), radiation (z), neutrinos {v), and a cosmological constant (A) as dark energy. 
For the basic background equations, we have from the energy-momentum tensor T for a 
scalar field, the scalar energy density Tq and the scalar pressure T^* are given by 

n = -p*o = -(^*o+^), (1) 

77 = P*o= (^*o-^)^i, (2) 

where the dots stand for the derivative with respect to the cosmological time and (5^- 
is the Kronecker delta. Thus, the cosmological Equation of State for the scalar field is 
P<jo = a;$oP$o with 

In order to solve the Friedmann equations with analytic methods with the approxima- 
tion m » H they performed a transformation and compared their result with numerical 
ones. Here the scalar field and the variables of the background depend only on time, e.g.. 

They computed the growth of the SFDM overdensities 6p^ in the linear regime, in 
this regime, the density contrast 6 = (5p$/p$o was much smaller than unity. It is believed 
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that the Universe was almost uniform after inflation, with a very small density contrast. 
As the Universe expanded, the small overdensities grew until they began to collapse, 
leading to the formation of structure in the Universe. Thus, only small deviations in the 
FLRW model are considered, so that they can be treated by linear perturbation theory. 
After introducing the perturbed metric tensor in the FLRW background, only scalar 
perturbations considered. In their paper they gave the equation of energy-momentum 
conservation and the Einstein field equations for the perturbed metric. Within the linear 
theory of scalar perturbations the evolution of the density contrast can be written as 

S + 3H (< > - < >] S ^ 2>(t)k < F<j, > - < Co > 



5p. 



where 



F$ = 1 + 

being (pk the gravitational potential. This equation differs from the density contrast 
equation for ACDM. However, in [76] they show that the extra terms F$ and G$ tend to 
the values of the standard equation of ACDM. Therefore the scalar perturbations in this 
model grow up exactly as in the ACDM paradigm. 

The evolution of the scalar perturbations in the nonlinear regime when 5 » \ was 
also studied. Here, an analysis was made within the framework of the spherical collapse 
model |91j . This formalism is very useful to understand the structure formation process 
in the Universe in the nonlinear regime. Focused on the era where the radiation density 
is equal to the SFDM density, the T « Tc, and therefore, it is expected that the scalar 
potential reaches the $q profile. They also studied if $^-dark matter perturbations (once 
the breaking of symmetry was achieved and the SF had reached its minimum ) where 
able to form bound structures as in the standard model. 

Following the same path, [114] obtained that for the matter dominated era the low-fc 
modes grow. When CDM decouples from radiation in a time just before recombination 
it grows in a milder way than it does in the matter dominated era (Fig. |3]). 
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ure 2: Evolution of the density contrast 6 for a perturbation with wavelength Afc ~ 2Mpc 



Although in general a scalar field is not a fluid, it can be treated as if it behaved like 
one and the evolution of its density can be the appropriate for the purpose of structure 
formation, because locations with a high density of dark matter can support the formation 
of galactic structure. 

In |114) they assumed that there was only one component to the mass density, and 
that this component was given by the scalar field dark matter. In this case the equation 
for the perturbations reads 



0, 



valid for all sub-horizon sized perturbations in the non-relativistic regime. 

It was shown that the scalar field with an ultralight mass of 10~^^eV simulates the 
behavior of CDM in a Universe dominated by matter when A = 0, because in general in a 
matter dominated Universe for low-fc, Vq (called the quantum velocity) tends to be a very 
small quantity tending to zero, so from the equation of the density contrast we could see 
that on this era we have the ACDM profile given by 



d^ 
di^ 



2H- 



J5 
di 



+ c 



AnGpo ]S = 



i.e., the SFDM density contrast profile is very similar to that of the ACDM model. Fig. [3j 
On the contrary for A 7^ both models have different behavior as can be seen from Fig. 
|4j results which show that linear fluctuations on the SFDM can grow, even at early times 
when the large-scale modes (small k) have entered the horizon just after aeq ~ 10~^, when 
it has decoupled from radiation, so the amplitudes of the density contrast start to grow 
faster than those for CDM around a ~ 10~^. Here an important point is that although 
CDM can grow it does so in a hierarchical way, while from Fig. [4] we can see that SFDM 
can have bigger fluctuations just before the ACDM model does, i.e., it might be that no 
hierarchical model of structure formation is needed for SFDM, and it is expected that for 
the non-linear fluctuations the behavior will be quite the same as soon as the scalar field 
condensates, which could be in a very early epoch when the energy of the Universe was 
about ~ TeV. These facts can be the crucial difference between both models. 
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Figure 3: Evolution of the perturbations for the CDM model (red dots) and SFDM model 
(green line) for k = 1 x 10^'^hMpc~^ . Notice how after the epoch of equality {a^q ~ 10^^) 
the evolution of both perturbations is identical, a = 1 today. In this case the self-interacting 
parameter is A = |114j 

Additionally, in '60j the growth of cosmological perturbations to the energy density 
of dark matter during matter domination was considered when dark matter is a scalar 
field that has undergone Bose-Einstein condensation. In this case, the inhomogeneities 
where considered within the framework of both Newtonian gravity, where the calcula- 
tion and results resulted more transparent, and General Relativity. The direction taken 
was again in deriving analytical expressions, which where obtained in the small pressure 
limit. Throughout their work the results where again compared to those of the standard 
cosmology, where dark matter is assumed pressureless, using analytical expressions to 
showcase precise differences. They also find, compared to the standard cosmology, that 
Bose-Einstein condensate dark matter leads to a scale factor, gravitational potential and 
density contrast that again increases at a faster rate. 

4.3 Galaxies and BEG dark matter mass constraints 

Magana et al. ,t3j considered a model where ultra-light bosons are the main components 
of the dark halos of galaxies. The main goal of this work was to constrain the mass of 
the scalar particles. They constructed stable equilibrium configurations of SFDM in the 
Newtonian limit to model the DM halo in UMi. They studied two relevant cases of SFDM 
halos: with and without self-interaction. 

Since galactic halos are well described as Newtonian systems, the work was done 
within the Newtonian limit. In this limit, The Einstein-Klein-Gordon equations for a 
complex scalar field $ minimally coupled to gravity and endowed with a SF potential 
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Figure 4: Evolution of the perturbations for the CDM model (red dots) and SFDM model 
(green line) for A; = 1 x W~^hMpc~^ and A 7^ and negative. Notice how after the epoch of 
equality {aeQ ~ 10^^) the evolution of both perturbations is now different from the one in 
Fig. [3j a = 1 today. In this case we can clearly see that the SFDM fluctuations grow faster 
than those for the CDM model. 



y($) = mj,$'^/2 + A$*/4, can be simplified to the Schrodinger-Poisson equations: 



ihdtip 



2m, 



A 



2m0 



(4) 
(5) 



where is the mass of the boson associated with the scalar field, U is the gravitational 
potential produced by the DM density core, A is the self-interacting coupling constant, 
and the field ip is related to the relativistic field $ through 



$ = e 



't/h 



UMi is a diffuse dSph galaxy located at a distance of 69 ± 4kpc from the Milky Way 
center and has a luminosity of Ly = 3 x IO^Xq. Its stellar population is very old with an 
age of 10-12Gyr. Dynamical studies suggest that UMi is a galaxy dominated by DM, with 
a mass-to-light ratio larger than 60Mq / Lq . Among the most puzzling observed properties 
of UMi is that it hosts a stellar clump, which is believed to be a dynamical fossil that 
survived because the underlying DM gravitational potential is close to harmonic. This 
condition is accomplished if the DM halo has a large core. 

In that work, it is mentioned that the most remarkable feature in UMi structure is 
the double off-ccntercd density peak. The second peak or clump is located on the north- 
eastern side of the major axis of UMi at a distance of ~0.4kpc from UMis center. The 
velocity distribution of the stars contained in the clump is well fitted by two Gaussians, 
one representing the background. The most appealing interpretation is that UMis clump 
is a disrupted cluster with an orbit in the plane of the sky, which has survived in phase- 
space because the underlying gravitational potential is harmonic, implying that the dark 
halo in UMi has a large core [5] 
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Figure 5: Snapshots of the clump in UMi galaxy, at t = 0, 5, and lOGyr. The clump is set on 
a circular orbit in the (x,y)-plane at a distance of = 0.39kpc from UMis center. The mass of 
the boson is = 10~^^eV and A = 0. The total mass of the galaxy is M = 9.7 x lO^M©. 



The fact that SFDM halos have cores might solve other apparent problems in dwarf 
galaxies. In their work the timing problem of the orbit decay of GCs (Globular clusters) 
in dwarf elliptical galaxies and dSph galaxies was considered. In fact, in a cuspy halo, 
GCs in these galaxies would have suffered a rapid orbital decay to the center due to 
dynamical friction in one Hubble time, forming a nucleated dwarf galaxy. For instance, 
under the assumption that mass follows light or assuming a NFW profile, Fornax GCs 3 
and 4, which are at distances to the center <0.6kpc, should have decayed to the center 
of Fornax in ~ 0.5-lGyr; this clearly represents a timing problem. Assuming a cuspy 
NFW halo, GCs 1, 2, 3 and 5 can remain in orbit as long as their starting distances 
from Fornax center are >1.6kpc, whereas one of them needs an initial distance >1.2kpc. 
However, there is no statistical evidence to suggest that the initial distribution of GCs 
is so different to the stellar background distribution. In addition, studies of the radial 
distribution of GCs in giant elliptical galaxies show that the distribution of metal-rich 
GCs matches the galaxy light distribution. Assuming that GCs formed along with the 
bulk of the field star population in dwarf galaxies, the probability that Fornax GCs were 
formed aU beyond 1.2kpc is ~ (0.03)^ ~ 2.5 x 10~^. Therefore, it is very unlikely that 
all the GCs in Fornax were formed at such large distances and even if they did, there is 
still a timing problem. 

The persistence of cold substructures in UMi places upper limits on m^. Using N- 
body simulations, it was found that the survival of cold substructures in UMi was only 
possible if 7710 < 3 X 10~^^eV in the A = case. On the other hand, by imposing a 
plausible upper limit on M, lower limits on where placed. All together, it was found 
that for A = 0, should be in the window 

0.3 X 10-^^eV < < 3 X IQ-^^eV. 

Since the timing problem of the orbital decay of the GCs in Fornax can be alleviated if 
< 1 X 10~^^ eV for A = 0, the most favored value resulted around (0.3—1) x 10~^^eV. 
For SFDM models with self-interaction, the upper limit on increases with A. 
Bosons of mass ^ 6 x 10~^^eV could account for the observed internal dynamics of UMi. 
In the limit A >> 1, it was found that m^/X < 0.55 x lO^eV* would explain the longevity 
of UMis clump and the surviving problem of GCs in Fornax. 

The window of permitted values for resulted quite narrow. Even so, it is remark- 
able that the preferred range for the mass of the boson derived from the dynamics of 
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dSph galaxies resulted compatible with those given by other authors to ameliorate the 
problem of overabundance of substructure and is also consistent with the CMB radiation 

[551[7ilM]. 

In a recent posting, Slepian and Goodman constrained the mass of bosonic DM us- 
ing rotation curves of galaxies, and Bullet Cluster measurements of the scattering cross 
section of self-interacting DM under the assumption that these systems are in thermo- 
dynamic equilibrium. If their assumptions are verified, repulsive bosonic DM will be 
excluded and, thereby, the only remaining window open is non-interacting bosons. Nev- 
ertheless, the static diffusive equilibrium between Bose-Einstein condensate and its non- 
condensated envelope, as well as finite temperature effects need to be reconsidered. In 
addition, other authors |103l I107| argue that scattering cross sections for bosonic DM are 
much smaller than those derived from the condition of thermodynamic equilibrium by 
Slepian and Goodman. 



4.4 SFDM density profiles and LSB galaxies 

An analysis of the Newtonian regime at Temperature zero can be found in Bohmer 
and Harko (2007), they assumed the Thomas- Fermi approximation which neglects the 
anisotropic pressure terms that are relevant only in the boundary of the condensate, the 
system of equations describing the static BEG in a gravitational potential V is given by 



= -P^V, (6) 



VV = 47rGp, (7) 

with equation of state 

p{p) = UoP^, 

where Uq = ^'^^ ° , p is the mass density of the static BEG configuration and p is the 
pressure at zero temperature, p is not the usual thermal pressure but instead it is produced 
by the strong repulsive interaction between the ground state bosons. 

In |104j , Robles and Matos found that the BEG dark matter model can give a density 
contrast profile consistent with RG's of dark matter dominated galaxies. The profile 
resulted as good as one of the most frequently used empirical core profiles, the pseudo 
Isothermal profile (PI) , but with the advantage of coming from a solid theoretical frame. 
In [104] The data was fitted within Ikpc and a logarithmic slope a = —0.27 ± 0.18 was 
found in perfect agreement with a core. They emphasized that the cusp in the central 
regions is not a prediction that comes from first principles in the GDM model, it is a 
property that is derived by fitting simulations that use only DM. For a detail discussion 
of the cusp and core problem in the standard model see [?] and references there in. They 
also explained an ambiguity in the usual interpretation of the core radius, they proposed 
a new definition for the core and core radius that takes away this ambiguity and that 
has a clear meaning that allows for a definite distinction of when a density profile is core 
or cusp. Using their definition they found the core radius in the BEG profile to be in 
most cases over 2kpc bigger than the core radius in the PI profile. They assumed a great 
number particles were in the ground state in the form of a condensate. This led to good 
results for their sample of galaxies, but it proved necessary to consider more than these 
simple hypotheses when dealing with large galaxies. 

The solution to the system above is 

PB{r)^pl^^ (8) 
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it can be seen that the BEC model satisfies p ^ r'^ near the origin, but a priori this 
does not imply consistency with observed RC's. Therefore, the profiles were fitted to 
thirteen high resolution RC data of a sample of LSB galaxies. The RC's were taken from 
a subsample of de Blok et al. (2001), galaxies that have at least 3 values within ^Ikpc 
where chosen, not presenting bulbs and the quality in the RC in Ha is good as defined in 
McGaugh S. S. et al. (2001). Because the DM is the dominant mass component for these 
galaxies they adopt the minimum disk hypothesis which neglects baryon contribution to 
the observed RC. In order to show that in LSB and dwarf galaxies neglecting the effect 
of baryons was a good hypothesis, they included two representative examples, see Fig. [6j 




Figure 6: Contribution of the baryons to the rotation curve for F568-3 and F583-1. We 
denote observed data by black dots with error bars, dark matter with blue asterisks, the disk 
with cyan squares and the the gas with magenta squared boxes. The figures on the left were 
fitted asumming the minimum disk hypothesis while the ones on the right are only the dark 
matter. In the fits shown are, BEC in solid line(red in the online version), PI dashed (green 
in the online version) and NFW double-dashed (black in the online version) profiles. 

For these galaxies the contribution of the gas was plotted, disc and dark matter 
separately. They did the fitting first considering the total contribution and then using 
only DM. They found no substantial difference in their values. 

As a second result and direct consequence of the core definition, they where able 
to obtain the constant value of /xq which is proportional to the central surface density. 
This result is one of the conflicts of the current standard cosmological model due to the 
hierarchical formation of galaxies. 

As the density profile eq. Q is not enogh to describe the large galaxies as discussed in 
|105j , Robles and Matos thus gave a physically motivated extension to the SFDM model 
that includes the DM temperature corrections to the first loop in perturbations. Their 
idea is to use the Z2 spontaneous symmetry break of a real scalar field as a new mechanism 
in which the early DM halos form. As stated earlier, when the real scalar field rolls down 
to the minimum of the potential, the perturbations of the field can form and grow. They 
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gave an exact analytic solution for an static spherically symmetric SF configuration, 
which in the SFDM model represents a DM halo. Their solution naturally presents a 
flat central density profilejust as eq.(|8]), but now it can accommodate more than just the 
ground state as the temperature T 7^ 0, in this way they solved previous discrepancies in 
rotation curve fits at T = 0, for instance, having a constant halo radius for all galaxies 
and the incapability to fit at the same time the inner and outermost regions of RC in 
large galaxies. Both issues were solved using this scenario which includes temperature of 
the DM and the exited states of the SF. 

The perturbed system of a scalar field with a quartic repulsive interaction but with 
temperature zero has been studied before [351 HH] ■ The study for the evolution of the SF 
with the temperature correction in a FRW universe is analogous. The metric tensor was 
written as g = Qq + Sg, where as always is the unperturbed FRW background metric 
and Sg the perturbation. The perturbed line element in conformal time 77 given by 

ds^ ^ a{rif[-{l + 2V')d?7^ + 2B,i drjdx' + (1 - 20)% + 2E„j dx'dx'] 

with a the scale factor, ip the lapse function, (j) gravitational potential, B the shift, and E 
the anisotropic potential. The energy-momcmtum tensor and the field where separated as 
T — To + 6T and ^{x^^) — $o('7) +(5$(a;^) respectively. As the linear regime was studied 
S^{x^) << $0(^7)1 the approximation ~ V{^q) could be made. They worked in 

the Newtonian gauge where the metric tensor g becomes diagonal and as a result, in the 
trace of the Einsteins equations the scalar potentials ip and are identical, therefore, ^ 
relates to the gravitational potential. 

The work was mainly focused in the galactic scale DM halos after their formation. 
Robles and Matos constrained themselves to solve the Newtonian limit of the perturbed 
KG equation, that is, 

□(5$ + ^[kl{T^ -Tl) + 12$2]J$ _ 4$o0 + ^[fc|(T2 - T^) + 4$2]$o0 = 

valid when $ is near the minimum of the potential and after the SB, where the SF is 
expected to be stable. Here Tc is the temperature of the symmetry break. 

Additionally to solving these two disagreements they mentioned why it does not seem 
necessary to include high amounts of feedback to fit and reproduce the inner core and 
wiggles found in high-resolution RC's,see Fig. [7] Also, this model can be tested with 
high redshift observations, the SFDM model predicts initial core profiles as opposed to 
the initially cuspy ones found in CDM simulatios which are expected to flatten due to 
redistribution of DM by astrophysical processes. 

Finally one conclusion is that if observations of more galaxies with core behavior are 
confirmed, this model can be a good alternative to ACDM. 

4.5 BEC dark matter and the power spectrum 

From another point of view in |120j it was again assumed that dark matter is composed of 
scalar particles that are able to form a Bose-Einstein condensate at some critical redshift 
Zcri but in this case it was used to study the matter power spectrum. 

After the BEC forms its effective pressure can assume a polytropic equation of state 
such as Phe ~ if an arbitrary non-linearity term is assumed, in this case 7. The exact 
value of 7 is defined by the non-linear contribution of the Gross-Pitaevskii equation which 
in its standard form leads to [18] 

Pbe — ^ Pbe- 
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Figure 7: Rotation curve of NGC 1560 and NGC1003 from [105j, where non zero temperature 
is considered. The wiggles are visible in the outer region of the RCs.The solid hne is the 
temperature corrected velocity profile and the dashed line is the Einasto fit |105] . 



In this case, the scattering length la and the mass of the dark matter particle again 
determine the dynamics of the fluid. Assuming that the condensate does not interact 
with any other form of energy, the above pressure, via the conservation balance, leads to 

^ Po 

where 

1.266 X r^beo X (m^/lmeV)-''' x (^^/lO^fm) 
^" ^ 1 + 1.266 X llbeo X (TOx/lmcV)-3 x (/„/109fm)' 

The current value of the scale factor a was taken as ao 1 and the current fractional 
density of the BEC dark matter as ^bea = Pfceo / Pc where pc is the critical density. When 
the above relations are combined the equation of state parameter of the BEC dark matter 
is obtained 

Po 

Wbe = —n . 

a* - Po 

If the inertial effects of the pressure become relevant as for example during the radia- 
tion phase or at the onset of the accelerated expansion, Newtonian cosmology fails and a 
more appropriate set of equations is needed. The inclusion of pressure in the Newtonian 
cosmology gives rise to the neo-Newtonian cosmology. In this case, the matter power 
spectrum is defined as always 

P(fc) = |,5b(z = 0;fc)|2, 
where (5;, (A:) is the baryonic density contrast calculated from equations 
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at the present time |120j . The baryonic agglomeration (5^, is supposed to be driven by the 
gravitational field which is sourced by all the forms of energy. In order to solve this set 
of equations it was needed to set the initial conditions for 5b and She and their derivatives 
at Zcr where the condensation took place. S {la,rnx) for each chosen couple 

of values (la,Wx) different initial conditions where needed. 

In [120] it was shown that if such phase transition occurred in the recent Universe 
this process would be able to leave small, but perceptible, imprints on the large scale 
structure perceptible in the matter power spectrum. Assuming la = lO^fm the BEC 
dark matter model does shows differences of the order of a few percents for masses 15-35 
meV. Adopting la — 10^°fm corrections of the same order where obtained for masses 
300-700meV. 

Although the BEC phase is shown to have a small influence on the matter power 
spectrum, a more quantitative analysis could be performed to estimate the preferred 
values of the model parameters. 

For the relevant parameter values studied in that work, the transition to the BEC 
phase was shown to occur at low redshifts. Since the standard cosmology remains un- 
changed before the CMB physics at the last scattering surface remained the same. 
However, the BEC dark matter would modify the gravitational potential just after Zcr 
while the speed of sound is nonzero leading to a contribution to the integrated Sachs- Wolfe 
effect. 



She 

(9) 



4.6 Vortices in BEC dark matter 

The conditions to form vortices in a SFDM/BEC halo have also been studied in |1031I125] . 
As it has been pointed out, these halos can be described as fluids, obeying quantum- 
mechanical fluid equations, so that the effects that make up this form of dark matter 
behave differently from standard ACDM, resulting in new effects with potentially observ- 
able consequences. The idea is that ultralight particles with m « leV will have very 
large de Broglie wave lengths which means that quantum statistical effects are important 
and macroscopic coherent lumps of matter can emerge. These light Bose particles will 
have a transition temperature to the condensed state that is of order ~ 2 • lO^K, which 
is the expected temperature in the Universe after about 1 second. 

There are essentially two limiting cases that can be considered. First, for quantum- 
coherence to be relevant on the scale of a halo of radius i?, the particle de-Broglie wave- 
length 



should be considered to be of the order of the halo size, XdeB S or ^Ise require 
^deB << R but with a strong repulsive self-interaction to hold the halo up against gravity. 
In the first case, if ~ Vmriai for the halo, this translates into a condition for the dark 



21 



matter particle mass, m>mH = 1.066 • IQ-^^iR/lOOkpcy^^^iM/lO^^MQy^^^eV cm^ 
where uih is a mass that depends on the properties of the halo and gn is determined 
by the density and the radius of the halo. In the second case, for the repulsive self- 
interaction pressure force to exceed the quantum pressure, it is required that g >> gn = 
2.252 • 10^'^4(i?/100kpc)(M/10i2Mo)-ieV cm^ If R is taken to be the radius of the 
virialized object supported against gravity by the dominant repulsive self-interaction, 
this imposes a condition on the particle mass given by m > ^Jlhg/gH■ 

However, it seems that rotating BEC haloes add new phenomenology, and the pos- 
sibility to distinguish this form of dark matter from other candidates. To this aim, in 
[103] they have studied the question of whether an angular velocity could be sufficient to 
create vortices in BEC/CDM cosmologies. As quantum fluid systems, BEC haloes can 
be modeled as uniformly rotating ellipsoids, with and without internal motions super- 
posed. To this aim, in |103j the authors derived equations which relate the eccentricities 
of haloes to their A-spin parameter. Once the latter is fixed, the eccentricities can be 
uniquely determined. They analytically studied necessary and sufficient conditions for 
vortex formation. In their results they found that vortex formation requires as a neces- 
sary condition that the halo angular momentum satisfies L > Lqm = NH^ which implies 
a lower bound on m/niH, i.e. on the dark matter particle mass. However, a sufficient 
condition for vortex formation could be established by an energy analysis, which aimed 
to find the conditions of when a vortex becomes energetically favored. 

They studied two classes of models for rotating halos in order to analyze stability with 
respect to vortex formation in two limits, one for L/Lqm >> 1 and for L/Lqm = 1, 
respectively. In what they called Halo-Model A {L/Lqm >> 1) these where modeled as 
homogeneous Maclaurin spheroids. The minimum angular momenta for vortex formation 
in this case was {L/LQM)crit = (5.65,4.53,4.02) for A = (0.01,0.05,0.1), respectively, 
which corresponded to a constraint on the particle mass m/mu > {m/mH)crit, where 
{m/mH)crit = (309.41,49.52,21.73), respectively. As long as m/mH satisfied this condi- 
tion, the strength of the self-interaction also satisfied the condition g/gn > {9/gH)crit, 
where {g/gH)crit = (1.02 • 10^,2549.24,454.54) for the same A-values, respectively. 

For Halo-Model B [L/Lqm — 1), which was an (n = l)-polytropic Riemann-S el- 
lipsoid, strictly irrotational prior to vortex formation, even L/Lqm = 1 resulted to be 
sufficient for vortex formation if the self-interaction strength was large enough. The 
condition L/Lqm = 1 fixed the value of m/mH for each A according to 



L m Kn 



-■1^1 



Lqm mn 10 (2 - e2)(l ~ elf/^{l - e2)i/3 

rr, ur /9R \^/^ / ^4 A 

m Kn I /-D12 \ /f, I ^1 \ ^1 



mHl0''[q,,) [^^Ail~el)) (l-e2)5/6(l-ei)V3' (10) 
and the condition of virial equilibrium 

y{x) = ^5(61,62)^^/^3;, 

thereby also fixing g/gn- For A — (0.01, 0.05, 0.1), these values where given by m/mn = 
(44.58,9.49,5.01) and g/gn = (1595.07,68.00,17.20), respectively Halo-Model B then 
made vortex formation energetically favorable for these values of m/mu and g/gn- They 
interpreted this to mean that, for L/Lqm > 1, vortex formation will also be favored, as 
long as g/gH > {9/ 9h)- Furthermore, any values of m/m^ and g/gu which satisfy the 
condition for vortex formation in Halo-Model A would automatically satisfy that found 
by Halo-Model B, which resulted less stringent but more accurate. 
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In conclusion they imagined vortex formation in BEC haloes composed of repulsively 
interacting particles as follows: If the angular momentum of a rotating BEC halo fulfills 
L < Lqm, no vortex would form, and the halo can be modeled by a mildly compressible, 
irrotational Riemann-S ellipsoid, which has a polytropic index of n = 1. For L = Lqm, 
the irrotational Riemann-S ellipsoidal halo can make a transition to a non-rotating, spher- 
ical halo with a vortex at the center if the self-interaction is strong enough. For a range 
of angular momenta fulfilling Lqm < L < 2Lqm, a central vortex can be expected but 
now with the excess angular momentum deforming the halo such that again a Riemann-S 
ellipsoid forms. Finally, if i >> Lqm, oblate haloes described as Maclaurin spheroids 
had a central vortex if m/mH > {m/mH)crit and g/gH ^ {g / gH)crit with the critical 
values given by Halo-Model A. Those critical values determined when a single vortex was 
energetically favored, but since L/Lqm » 1, it is also possible that multiple vortices 
could form. 

From another point of view, in [125] it was assumed that the particles where non- 
interacting and therefore only gravity acted on the system. Following [5S], the authors 
based on Jeans instability analysis to estimate their parameters. The growing mode under 
gravity was given by e"*"* with 7^ = AnGp, whereas the free field was supposed oscillatory; 
g-iEt ^[^Yi E = k'^ /2m. The latter was then written as e'^* with 7^ = — (fc^/2rn)^. 
Noting that this is like normal Jeans analysis with sound speed — k/2m then 7^ = 
AirGp — (fc^/2m)^. Setting this to zero, for the Jeans scale they found 

rj = 2^/kj = 7r3/4(Gp)i/4m-i/2 = h-,m-^'\p/ p.y^l^a^h'r^'^^c, 

where 77122 — Tn/lO^^^eV, and the background density is pb — 2.8 • lO^^r^m/i^Ai^QMpc"^. 
It is supposed that below the Jeans scale the perturbations are stable and above it they 
behave as ordinary CDM f55|. As always, the stability below the Jeans scale was guar- 
anteed by the uncertainty principle. If the particles are confined further, their momenta 
increases and oppose gravitational contraction. 

Moreover, in (1251 the author also considered the suggestion that superfluid BEC dark 
matter in rotation would likely lead to vortices as seen in atomic BEC experiments. As 
already noted in [40], BEC dark matter with self-interactions could actually constitute 
a superfluid. In the case of a repulsive self-interaction in [125j the author too argued 
that the vortex size should be determined locally by the coherence length. This means 
that there would be two scales in the problem: A galactic one, given by the de Broglie 
wavelength from the tiny mass, and a sub galactic one determined by the mass and the 
two-body interaction strength. He also explored the consequences of self-interactions 
on the virialization of gravitationally bound structures and found almost no effects for 
reasonable values of 771 and a. Here the case considered was also that of a quartic self- 
interaction Ac/)"*. 

In [IH] there was a brief discussion of the effect on the Lane-Emden equation, whereas 
in other works a BEC of axions with a single vortex arising from global rotation in the 
early Universe has been considered. This latter scenario is, however, less likely to occur 
since the global rotation rate of the Universe can be estimated from various observations 
and is very small but nevertheless non-zero. 

Under the assumption of dark matter being an ultra-light BEC, the rotation of spiral 
galaxies would cause vortex lattices to form. In |125j . the author also considered possible 
effects of sub galactic vortices in the dark matter on the rotation velocity curves of viri- 
alized galaxies with standard dark matter halo profiles. He found that one can actually 
get substructure in the rotation curves that resemble some observations, but that this 
required large vortex core size and small vortex-vortex distances. The mass and interac- 
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tion strength needed to realize this where found to be fine-tuned, but could possibly be 
accommodated in more general setups. 

If dark matter contains a component of condensed BEC particles that is superfluid and 
if the halos are rotating then it is not inconceivable that there can be vortex formation. 
However, the quantized vortex discussion makes an important assumption about the 
coherence length, ^, entering fie (angular velocity of the halo) in 

h , /1.46i?\ 6.21 -10-1^ /l.46i?\ 

Here ^ is taken to be of kpc size. It then results that with no self-interaction there is only 
the gravitational scale fi^/GMm? available, which becomes of galactic size for masses 
m ^ 10~^^eV. However, when including self-interactions through the scattering length a, 
there is also a scale given by ^ = l/VSTrart, which is the usual Gross-Pitaevskii coherence 
length. The latter coincided with the characteristic length over which the density is 
expect to go to zero in a vortex. 

If the additional assumption that the vacuum expectation value, (j)^), arises from a 
mechanism that preserves parity, the interaction term results 

This terms is of course merely the standard interaction term in the Gross-Pitaevskii theory 
of interacting condensed bosons. Therefore he also concludes that the self-interacting 
scenario emerges from this procedure. 

These investigations and simple numerical experiments pointed to an interesting effect 
that could arise from bosonic dark matter. However, to fully explore the influence that 
vortex lattice formation and stellar feedback on structure formation has in luminous 
matter, an ultralight BEC dark matter component in large N-body simulations should 
be considered. 



4.7 Bose dark matter additional constraint 

In [107j BE condensation inside the primeval fireball, at zero-order in perturbation theory 
has been studied. Here, the process of condensation was considered to be driven by self- 
interactions of high-energy bosonic particles. 

In this work it was found that in the instantaneous decoupling approximation, the 
subsequent evolution of the full bosonic system was only affected by the expansion rate 
of the Universe and small gravitational instabilities. 

The evolution of bosonic DM after decoupling was analyzed as follows: their velocity 
and temperature affected only by the expansion rate of the Universe. 

In 

^cPcr ^HPcr 



the temperature was needed in order to calculate the mass of the bosonic DM. Here 
fic and fin represent the content of bosonic CDM and bosonic HDM (hot dark matter) 
respectively, Uc is the number density of condensed bosons and ht the number density 
of thermal bosons. 

As an additional remark, they also address a bound on the strength of the bosons 
self- interactions . 
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In their study, the bosonic DM parameters where addressed as the mass, to, and the 
factor gx (amount of degrees of freedom); where bounds on their values have been obtained 
from a statistical analysis of cosmological data. The constraints found for the temperature 
of the boson gas T^^ = 2.14±0.02K, and for the boson-antiboson gas T^^ = 1.91±0.05K. 

Finally, from a similar analysis for fermionis they found bounds in the sum of neutrino 
masses and the number of extra relativistic species, ^ < 0.45eV, N ~ 1.10 ± 0.18, in 
concordance with some previous reports. 

In summary, they presented a generic study of DM based on BE condensation, from 
which, CDM and HDM, result intrinsically related. 



4.8 SFDM and black holes 



The rapid decay of the energy density of the scalar field for the case of super-massive 
black holes, indicates that scalar fields may not be maintained around a black hole during 
cosmological time scales in the whole space; so either the scalar field gets accreted or it 
escapes through future null infinity. The fact is that when a Schwarzschild black hole that 
is asymptotically flat is considered, scalar fields tend to vanish from the spatial domain. 

An appropriate coordinate system to study the propagation of scalar fields in the 
Schwarzschild space-time is using hyperboloidal slices, because it has been seen that 
such slices reach future null infinity instead of spatial infinity, which results as a natural 
boundary for a wave-like processes, including electromagnetic fields and gravitational 
radiation. 

In this case the Klein-Gordon equation for a scalar field (pr can be written as [30] : 



□02 

v = 



dV 



— 0, where Q^t- 



=^ dfj. W~99'^'^(^v4't] , with a potential of the form 



where has the units of mass. Then the KG equation in the 



conformal metric can be expressed as: 
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0, 



provided the relationship between the physical scalar field 0t and the conformal scalar 



field 4>T to be (px = 4>t/^- Here R 



12Q. 



[r + (2r — 1)] is the Ricci scalar of the conformal 



metric and □ = V^V^ corresponds to the conformal metric. The case for which this last 
equation results conformally invariant corresponds to the zero mass case tjib — 0. 

In their work jSOj, used initial scalar field profiles with m\ = 0.0,0.1,0.2. In order 
to explore the parameter space, various values of the amplitude A = 0.01,0.1,1.0 and 
different widths of a Gaussian pulse, which in physical units corresponded to cri — 0.5, 
(72 = 1 and tTs = 5 to the right from rg = 0.8. This range of parameters allowed the 
authors to cover length ranges that involve Compton wave-length related to effects of 
interaction and reflection that may involve reflection and absorption effects. 

In their work, [30 , one of the conclusions was that one potential ingredient that would 
help at maintaining massive scalar fleld densities during longer times is the rotation of 
the black hole and also of the scalar field. It would then be of major interest to use 
foliations that approach future null infinity and study if the same effects occured and also 
the study of scalar field configurations with non-zero angular momentum. In fact, the 
results in their work where obtained assuming the maximum cross section of accretion 
due to the spherical symmetry, which in turn worked as upper bounds accretion rates in 
more general cases. 

Yet another possibility results on considering solutions that asymptotically may con- 
tain a cosmological constant, which would be appropriate if a background energy density 
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in the universe where assumed. This would imply that black hole candidates should 
not be considered to be asymptotically flat. Other possibilities may include black hole 
candidates of a different nature like boson stars. 



4.9 Threats to scalar field dark matter, black holes? 



As mentioned in the last subsection, the existence of a long-lasting scalar field configura- 
tions surrounding a black hole have been studied [30| I13j . Another motivation for these 
studies is the possibility that super-massive black holes at galactic centers may represent 
a serious threat to the scalar field dark matter models. 

As a first step, a relatively simple model has been considered [13] ■ There, stationary 
scalar field configurations to the Klein-Gordon equation 

(□-^^2)0 = 0, 

where looked for on a Schwarzschild space-time background. 



ds^ = -N{r)dr + 



N{r) 



+ r^dQ.\ N{r) := 1 - 2M/r, 



with the d'Alambertian operator defined as □ := {^l \/—g)d^{^J—gg^^'^d^), M being the 
mass of the black hole and dil^ := dO^ + siv^Odip^ the standard solid angle element. So 
far, the case has been restricted to the case of a canonical, massive, non self-interacting 
minimally coupled scalar field </>. With these conventions (j) results dimensionless, while 
jjL has dimensions of length^ ^. The associated quantum mechanical "mass" of the scalar 
field given by h^. 

In order to look for the stationary solutions of 



dr ^ ^ dr 



-Uiifi,M;r) 



Iplm = 0, 



a further decomposition of the functions ipimit, was done into oscillating modes of the 
form: 

with uiim a real frequency and uimir) a complex function of r in the interval (2M, oo). 
Although stationary solutions where found for the scalar field. 



HI' 



N{r) 



2M 



du 
dr 



-Uli^l,M■,r)\u\^ dr 



they have been shown to be unphysical, in the sense that their energy density integrates 
to infinity in a compact region just outside the event horizon. However, there seem to 
exist long-lasting, quasi-stationary solutions of finite energy, which are found by evolving 
initial data that was constructed by slightly modifying a particular subset of the stationary 
solutions. The solutions found so far show as an overall behavior an exponential energy 
decay, caused by scalar field leaking into the black hole, that in some cases can be very 
slow. 

The stationary solutions where obtained by solving a time-independent Schrodinger- 
like equation 



dr* 



Veffin 



u{r*) — uP'u{r*), 



oo < r < oo. 
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with an effective potential 



Veffir*) ■.^N{r)Ui{fi,M;r), r = r{r*). 

and fience tfiey wfiere cfiaracterized only by the properties of the potential. This fact is 
strictly true for the stationary solutions, but interestingly it was found that the quasi- 
stationary solutions, for which the Schrodinger-like equation no longer holds, can also 
be characterized by the properties of that same effective potential. Then, the cases of 
interest where those in which the effective potential contained a local minimum given by 
the condition 

M^i^r^ - l{l + l)r^ + 3M{P + I ~ l)r + 8M^ = 0. 

The existence of this minimum then depended solely on the combination of the parameters 
and I. Although none of the possible forms of the effective potential proposed allowed 
for bound states, the existence of the well was enough to allow for resonant states, which 
where the ones that where useful in constructing initial data that could give rise to long- 
lasting quasi-stationary configurations of finite energy. 

As mentioned by the authors, it could be objected that in order to obtain the men- 
tioned quasi-stationary solutions very particular initial data should be constructed. How- 
ever, it seems that the crucial factor is the existence of the potential well. Even when 
starting with modified stationary solutions that are not resonant, after an initial abrupt 
energy loss, the late time behavior observed resulted very similar to that of the resonant 
quasi-stationary solutions. These solutions seemed to evolve as a combination of the res- 
onant modes. As mentioned before, the value of /i for scalar field dark matter models is 
expected to be given approximately by = lO^^^eV in physical units, which gives rise 
to effective potentials with a local minimum for all values of I. 

When evaluating the characteristic times of the solutions it was found that the longest 
lasting configurations could last for thousands of years. Although cosmologically this is 
very short time, for technical reasons, the authors where only able to study cases with 
relatively large values of M^. Noting how fast the characteristic times seemed to increase 
with decreasing fi, it resulted unreasonable to expect that configurations with ~ 10~^ 
could last for cosmological time-scales. 

Some aspects of this study are still open for some improvement. For example, a self- 
gravitating scalar field has not yet been considered. Second, much smaller values of the 
parameter /x, and much larger scalar field distributions, could be needed to do a more 
realistic representation of dark matter halos. The main difficulty in dealing with such 
configurations result in handling the very different scales numerically. Third, besides 
studying possible quasi-stationary or long-lasting configurations with an already existing 
black hole, it would also be interesting to consider more dynamical scenarios such as the 
formation and/or growth of the black hole and the possibility of survival of the scalar 
field afterwards. The results presented by Barranco et al. in [13] seem to indicate that 
it is indeed possible for scalar field halos around super-massive black holes to survive for 
cosmological time-scales. 

5 Conclusions 

In this work we have revisited an alternative DM paradigm of the Universe known as 
scalar field dark matter or Bose-Einstein condensate dark matter model. In this model a 
fundamental scalar field plays the role of dark matter. 
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We have reviewed a large number of recent quantitative and qualitative results aimed 
at explaining a variety of trends seen in the SFDM/BEC model. These trends include 
a brief description of Bose-Einstein condensates as dark matter, analysis and growth of 
its cosmological perturbations, its effect on the dynamics of galaxies, density profiles and 
the mass power spectrum, among others. 

The key parameters of an ultralight EEC dark matter model are naturally shown to 
be the mass of the boson, which must be extremely small and, for the self-interacting 
scenario, the coupling strength of two-body repulsive interaction among the condensed 
particles. Dark matter is then suggested to arise from a single scalar field coupled to 
gravity undergoing a spontaneous symmetry break and hence rolls to a new minimum 
which gives a new vacuum expectation value. The breaking of symmetry can be done via 
a Ginzburg-Landau potential with quadratic and quartic terms. These give the mass and 
the interaction terms to the scalar field. 

In the cosmological regime, it has been shown, that the SFDM/BEC model with 
an ultralight mass of 10~^^eV mimics the behavior of the cosmological expansion rate 
predicted with the ACDM model. 

Another interesting cosmological behavior of the SF indicates that their scalar fluctua- 
tions can be appropriate for the purpose of structure formation, mainly because overdense 
regions of SFDM/BEC can produce the formation of galactic structure. Thus, the stan- 
dard and the SFDM/BEC models can be contrasted in their predictions concerning the 
formation of the first galaxies. If in the future we observe more and more well formed 
and massive galaxies at high redshifts, this could be also a new indication in favor of the 
SFDM/BEC paradigm. 

Much of the interest in this model has focussed on its ability to predict and agree 
with observations of the existence of dark matter and its capability to compete with the 
standard model of cosmology ACDM. Most of the themes wc have described attempt to 
provide a wider view of what the model is and what is its current status. There are many 
other works in this field of research that are not less important, however, we trust that 
the interested reader will be able to deepen its knoledge with the references cited in this 
article. 

Clearly there are a large number of competing models to describe the dark matter of 

the Universe, and slight differences such as vortex formation or small shifts in cosmolog- 
ical perturbations may not be enough to decide which one is the best, however, detail 
observations give us a means to discard several of them and identify the ones that can 
stay. 

Finally, the SFDM/BEC model can have important implications in the nature of 
dark matter in the Universe. Additional work is needed if we are to fully understand 
this model, it would be desirable to have a unify framework that involves rotation curves, 
vortex formation, if present, and black hole effects. With all these intriguing results, 
we consider that the SFDM/BEC should be taken as a serious alternative to the dark 
matter problem in the Universe. The observational evidence seems to be in favor of some 
kind of cold dark matter, if we continue to observe even more galaxies at higher redshifts 
and if higher resolution observations of nearby galaxies exhibit a core density profile, this 
model can be a good alternative to ACDM. We expect that future observations in galaxy 
surveys can get us closer to the nature of dark matter. 
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